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ABSTRACT 

Supermassive stars (SMS) are massive hydrogen objects, that slowly radiate their gravitational 
binding energy. Such hypothetical primordial objects may have been the seed of the massive 
black holes (BHs) observed at the centre of galaxies. Under the standard picture, these objects 
can be approximately described as n = 3 polytropes, and they are expected to shine extremely 
close to their Eddington luminosity. Once however, one considers the porosity induced by in- 
stabilities near the Eddington limit, which give rise to super-Eddington states, the standard 
picture should be modified. We study the structure, evolution and mass loss of these objects. 
We find the following. First, the evolution of SMSs is hastened due to their increased energy 
release. They accelerate continuum driven winds. If there is no rotational stabilization, these 
winds are insufficient to "evaporate" the objects, such that they can collapse to form super- 
massive BHs, however, they do prevent SMSs from emitting a copious amount of ionizing 
radiation. If the SMSs are rotationally stabilized, the winds "evaporate" the objects until a 
normal sub-Eddington star remains, having a mass of a few lOOA/©. 
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1 INTRODUCTION 

Super-massive stars (henceforth SMSs) are hydrogen objects with 
large masses, in the range of 10'* — Kf Mq . Their state is described 
as a hydrostatic equilibrium between their self gravity on one hand, 
and radiation pressure on the other. Since radiation leaks out as the 
objects shine, the thermal energy is replenished through the slow 
release of the gravitational binding energy. This is unlike the less 
massive stars, where the central temperatures and densities are high 
enough to ignite thermonuclear reactions, the energy of which sus- 
tains the hydrostatic equilibrium. Thus, SMSs look like massive 
radiatively supported stars, but without nuclear reactions (Wagoner 
1969; Weinberg 1972; Shapiro & Teukolsky 1983). 

Under the standard picture (summarized in §1.1), SMSs are 
expected to shine at nearly the Eddington limit. However, it is now 
known that objects shining near this luminosity should develop 
instabilities, reduce the effective opacity in their atmosphere and 
attain super-Eddington luminosities (Shaviv 2001a,b). SMSs are 
expected to be no different, and should be described with super- 
Eddington states. 

The main qualitative difference introduced by these super- 
Eddington states, is the acceleration of a continuum driven wind. 
This high mass loss rate may affect their evolution. It is also re- 
sponsible for a larger photosphere which may affect the amount of 
ionizing radiation emitted by the SMSs. These effects are studied 
in the present work. 

Two further complications should be considered. First, it was 
demonstrated by Fowler (1966) that a dynamically insignificant 
amount of rotation (namely, rotation for which the star is still close 
to being spherical) can stabilize the SMS against the GR collapse 



into a massive black hole. In the present analysis, we will consider 
both cases — with and without rotational stabilization. 

The second complication was proposed by Begelman (2010), 
who considered a more realistic model for the formation of SMSs. 
It was argued that under a natural scenario, SMSs can form with 
an entropy inversion, such that the outer layers which are accreted 
later, have a higher entropy. Such SMSs are not n = 3 polytropes 
anymore. Instead, they have a small convective core surrounded by 
a massive convectively stable envelope. We will not consider this 
type of model here. Just as the standard model is modified here to 
consider the super-Eddington states, Begelman's scenario should 
be likewise modified. We defer this to subsequent work. 



1.1 Super Massive Stars: The Standard Picture 

To first approximation, SMSs can be modeled as pure polytropes 
(e.g., Shapiro & Teukolsky 1983). That is, they can be approxi- 
mated using P = Kp^ . Since radiation pressure dominates the 
gas pressure in SMSs, the total pressure is approximately given by 
the radiation pressure. Together with equation of hydrostatic equi- 
librium, 



dP _ GM{r)p 
dr r^ 

and equation of radiative transfer 
2 c d 



L 



-^-'^^3^d^(«^') 



(1) 



(2) 



it is easy to see that the star radiates at its Eddington limit. 
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Once we will allow the opacity to decrease in the super-Eddington 
states, we see that the object's luminosity will necessarily increase. 



7.7.7 Stability of non-rotating SMSs 

One of the most interesting aspects of non-rotating SMSs, is their 
instability towards gravitational collapse, which arises from cor- 
rections due to General Relativity. We succinctly summarize the 
conditions under which this instability arises (for more details see 
Shapiro & Teukolsky 1983). These are important as they will be re- 
quired once we will describe the evolution of the SMSs. The New- 
tonian energy of a polytrope is given by 



E = / 47rr ( u + 

'0 



GMp \ 



dr. 



(4) 



where u is the internal energy density, u = P/{T — 1) and V = 
4/3 is the polytropic index. This integral is solved for an n = 3 
polytrope to get 



E = kiKp-J'^M - kiGpl'^M^'^, 



(5) 



where fci and ^2 are constants, and pc is the central density of the 
SMS. For a given mass, the equilibrium configuration occurs when 
dE/dpc ~ 0. This condition of equilibrium provides a relation 
between the total mass and the specific entropy of radiation. 



Srad — 



4m H a /SttG 



(6) 



Higher order coiTections to the energy can be obtained by 
adding the effects of gas pressure to the internal energy and 
the post-Newtonian (i.e., GR) corrections to the gravitational en- 
ergy. The total energy of the star including these corrections is 
(Shapiro & Teukolsky 1983, chap. 6) 



.1/3 



5/3 1/3 



E = AMpt,'-' - BAV'-'py'' + CMpy-'lnp 



.1/3; 



D^f'^pl'^.^!) 



where A, B, C and D are constants, (the terms proportional to C 
and D are derived from the additional effects of gas pressure and 
post-Newtonian corrections). The condition for equilibrium config- 
uration dE/dpc ~ 0, gives the central density as a function of the 
mass and the specific entropy. General relativity instabilities ap- 
pears when d^E/dpc = 0. Differentiating eq. (7) twice, yields: 



0: 



-CMp: 



-2/3 



DA'f^'p:'^^. 



(8) 



Thus, there is a critical value for the density, above which the star 
becomes unstable. This analysis will not change once we introduce 
the existence of super-Eddington states. 



7.7.2 Stabilization by rotation 

The General Relativistic instability disappears if the SMS is rotat- 
ing fast enough (Fowler 1966). In this case, the Newtonian term for 
the rotational energy should be added to the total energy (eq. 7). It 
is given by 



*= i / r'^uj'^sin{e)pdV, 



(9) 



where u is the angular velocity and 6 is the polar angle measured 
from the axis of rotation. For a constant^ angular momentum $, the 
critical value of $ which stabilize the SMS against GR instability 
is given by 



—-r ~ 3.6 X 10 -n—crn sec , 
M Mq 



(10) 



this value is relatively small even for a SMS with M = 10 M©. If 
$ > <l?cr, rotation prevents the gravitational collapse. 



2 BACKGROUND: SUPER-EDDINGTON STATES 

Before proceeding to construct super-Eddington SMSr models, we 
begin by reviewing the relevant physics pertaining the emergence 
of super-Eddington states. These include three particular elements. 
First, the rise of inhomogeneities was shown (Shaviv 1998) to re- 
duce the effective opacity. Second, once a super-Eddington state 
arises, strong continuum driven winds are accelerated. Last, if the 
wind mass loss is too large, wind stagnation and a photon-tired 
state arises (Owocki & Gayley 1997). In it, a layer is formed in 
which strong shocks mediate a high energy flux without an exces- 
sive mass flux. These components are necessary building blocks 
for the super-Eddington states, and we therefore review them be- 
low. Two more examples where this theory is applied can be found 
in Dotan & Shaviv (201 1) and Dotan et al. (201 1). 



2.1 The rise of super-Eddington states 

According to common wisdom, objects cannot shine beyond their 
classical Eddington limit, T^Edd, since no hydrostatic solution ex- 
ists. In other words, if objects do pass £Edd, they are highly dy- 
namic. They have no steady state, and a huge mass loss should 
occur since their atmospheres are then gravitationally unbound and 
they should therefore be expelled. Thus, astrophysical objects ac- 
cording to this picture, can pass £Edd but only for a short duration 
corresponding to the time it takes them to dynamically stabilize 
once super-Eddington conditions are forced. 

For example, this can be seen in detailed ID numerical simu- 
lations of nova thermonuclear run aways, where novae can be in a 
super-Eddington state but only for several dynamical time scales 
(e.g., Prialnik c& Kovetz 1992). However, once they do stabilize, 
they are expected and indeed do reach in the simulations a sub- 
Eddington state. Namely, we naively expect to find no steady state 
super-Eddington atmospheres. This, however, is not the case in na- 
ture, where nova eruptions are clearly super-Eddington for dura- 
tions which are orders of magnitude longer then their dynamical 
time scale (Shaviv 2001b). This is exemplified with another clear 
super-Eddington object — the great eruption of the massive star tj- 
Carinae, which was a few times super-Eddington for over 20 years 
(Shaviv 2000) 

The existence of a super-Eddington state can be naturally ex- 
plained, once we consider the following. 

(i) Atmospheres become unstable as they approach the Edding- 
ton limit. In addition to instabilities that operate under various 
special conditions (e.g.. Photon bubbles in strong magnetic fields, 
Arons 1992; Gammie 1998; Begelman 2002, or s-mode insta- 
bility under special opacity laws, Glatzel 1994; Papaloizou et al. 



1 For the more general case of non-constant angular momentum see Fowler 
(1966). 
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1997), two instabilities operate in Thomson scattering atmospheres, 
(Shaviv 2001a). It implies that all atmospheres will become mista- 
ble already before reaching the Eddington limit. 

(ii) The effective opacity for calculating the radiative force on 
an inhomogeneous atmosphere is not necessarily the microscopic 
opacity. Instead, it is given by 



off 



{F)y 



(11) 



where { ) y denotes volume averaging and F is the flux. The situa- 
tion is very similar to the Rosseland vs. Force opacity means used in 
non-gray atmospheres, where the inhomogeneities are in frequency 
space as opposed to real space. For the special case of Thomson 
scattering, the effective opacity is always reduced. 

Thus, we find that as atmospheres approach their classical Ed- 
dington limit, they will necessarily become inhomogeneous. These 
inhomogeneities will necessarily reduce their effective opacity such 
that the effective Eddington limit will not be surpassed even though 
the luminosity can be super-classical-Eddington. This takes place 
in the external part of luminous objects, where the radiation dif- 
fusion time scale is shorter than the dynamical time scale in the 
atmosphere. Further inside the atmosphere, convection is necessar- 
ily excited such that the total energy flux may be super-Eddington, 
but the radiative part of it is necessarily sub-Eddington with the 
convective flux carrying the excess (Joss et al. 1973). 



2.2 Super-Eddington Winds 

The atmosphere remains sub-Eddington while being classically 
super-Eddington only as long as the inhomogeneities comprising 
the atmosphere are optically thick. This condition will break at 
some point where the density is low enough. At this height, the 
effective opacity returns to its microscopic value and hence the ra- 
diative force becomes super-Eddington again. Above this point we 
obtain a thick wind. Because it is optically thick, the conditions at 
the wind affect the structure at its base. 

At the critical point, the radiative and gravitational forces 
balance each other. This point will coincide with a sonic surface 
(where the mass loss velocity equals the local speed of sound) for a 
steady state wind. This allows us to obtain the local mass loss rate 
which is given by 



rhw = 47rrsPcrit«s(r-s) = const. 



(12) 



where r^ is the critical (sonic) radius, pci-it is the density at this 
radius and Vs is the local isothermal speed of sound. 

Based on the fact that instabilities develop structure with a typ- 
ical size comparable to the density scale height in the atmosphere, it 
is possible to estimate the average density at the sonic point (Shaviv 
2001b). Using this density the mass loss can be estimate to be 



W- 



(13) 



where W is a dimensionless wind "function". In principle, >V can 
be calculated from first principles only after the nonlinear state 
of the inhomogeneities is understood. This however is still lack- 
ing as it requires elaborate 3D numerical simulations of the non- 
linear steady state. Nevertheless, it can be done in several phe- 
nomenological models which depend on geometrical parameters 
such as the average size of the inhomogeneities in units of the 
scale height (/3 = d/lp), the average ratio between the surface 
area and volume of the blobs in units of the blob size (H), and 
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Figure 1. The general structure of a SMS above the Eddington limit. In 
regions where the density is high enough, the Eddington luminosity cannot 
be surpassed. Deep inside the star, convection will necessarily arise as the 
Eddington luminosity is approached, keeping the luminosity below £Edd- 
Further out, it was discovered that instabilities will arise as L approaches 
^Edd (£g- Shaviv 2001a), and that the ensuing nonlinear structure will 
have a reduced effective opacity (Shaviv 1998). However, as the overall den- 
sity declines with height, the clumps become optically thin and so are ex- 
posed to the full radiative flux, thereby increasing the net force above grav- 
ity and introducing a net supersonic outflow, forming a continuum-driven 
stellar wind. 



the volume filling factor q of the dense blobs. For example, in 
the limit in which the blobs are optically thick, one can show 
that W = 3H/3207a/3(l - af (Shaviv 2001b), with v being 
the ratio between the effective speed of sound in the atmosphere 
to the adiabatic one. Thus, W depends only on geometrical fac- 
tors. It does not depend explicitly on the Eddington parameter 
r = L/C-Edd- It yields typical values of >V ~ 1 - 10. 



3 THE MODEL 

We now proceed to describe the underlying physics necessary for 
the construction of a SED SMS. This can be divided into three 
parts. The first is a description of the internal structure of the SMS. 
Although it can be approximately described by a polytrope, we will 
do so with the more complete set of equations. The second part is 
that of the porous atmosphere which forms the outer static layer of 
the object. This layer is responsible for the reduced opacity which 
allows for the emergence of the super-Eddington state. Above this 
layer, we find the third part, which is that of a continuum driven 
wind. Because this wind is optically thick, it affects the structure, 
and later the appearance of the SMSs. 

The general structure is summarized in fig. I , and further elab- 
orate below. 



3.1 The Stellar Structure Equations 

3.1.1 Convective Core 

The convective zone begins at the center of the star and occupies 
the bulk of the hydrostatic radius and most of the SMS mass. In 
this region, the first two equations describing the stellar structure 
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(in radial coordinates, using standard notation), are the equation of 
hydrostatic equihbrium 

dP _ Gm{r)p 
dr r^ 

and the equation of mass conservation, 
dm 



(14) 



dr 



■ 4nr p. 



(15) 



The third is the equation of state, which is obtained as a combina- 
tion of an ideal ionized gas and black body radiation. 



S prrip 



(16) 



where p = {2X + 1.5Y)~^ is the mean molecular weight, and X 
& Y are the hydrogen and helium mass fractions respectively. 

The last equation describes energy transfer inside the star. 
Under conditions prevailing in SMSs, the relevant energy transfer 
takes place through convection. The temperature gradient is hence 
given by 

dT_ _ 7 - 1 T dP 
dr ^ 7 PdT' 



(17) 



where 7 is the adiabatic index. Convection is present if the standard 
Schwarzschild criterion is satisfied, but although this condition is 
satisfied all the way out to the stellar surface, at a certain radius the 
energy transport becomes purely radiative, as convection becomes 
inefficient. For convection to be efficient, the convective flux must 
be smaller than the maximum possible, which is given by 



iconv.max = 47rr pV ^ , 



where Vs is the adiabatic speed of sound. 



(18) 



(19) 



3.1.2 Atmosphere 

If the density is too low, convection becomes inefficient, the energy 
transfer takes place through radiation diffusion, and the radiative 
luminosity, Lrad becomes super-Eddington. The equation of radia- 
tive transfer is modified to take into account the porous nature of 
the atmosphere: 

dT _ ZKesLp 
dr 16nacr'^T^ ' 

where Kch is the effective opacity. As described in §2.1, when the 
flux approaches a critical value, the atmosphere develops inhomo- 
geneities and the gas becomes porous, such that the radiative force 
exerted on the gas is reduced. We assume that the relation between 
the effective Eddington factor Fefi = L/ Lea and the classical Ed- 
dington factor F = L/ C^dd is empirically given by 

'-^ 

F for F<Fcrit. 



refi 
refi 



for F > Fcrit, 



(20) 



Fcrit is the critical F above which inhomogeneities are excited, so 
the effective opacity for F > Fcrit can be written as 



Koff = KTh ( 1 - 7Tb ) /r- 



A 



(21) 



Since we expect a continuous Fefi, A, B and Fcrit satisfy the equa- 
tion Fcrit = 1 — ^/F^it. From theoretical considerations, we 
take Fcrit ~ 0.8 (Shaviv 2001a), though we have checked also 
Fcrit = 0.5, as elaborated in §5.4. This implies a relation between 
the normalization constant A and the power law B, that is given by 



^ = (1 - FcHt)Fc°it 



(22) 



The other equations describing the structure of this region is the 
equation of hydrostatic equilibrium, 

dP 

^ = -9P, (23) 

where g = GM/r^ . To a very good approximation, the mass in the 
atmosphere is negligible, such that we can safely assume that 

dM 



dr 



= 0. 



(24) 



3.1.3 Continuum Driven Winds 

As we described in §2.2, the average density decreases as we move 
up the atmosphere. At some point, the typical optical depth of the 
nonlinear structure would decrease below unity. Once the typical 
structure becomes transparent, radiation cannot be funneled around 
anymore, and the effective opacity returns to its microscopic value. 
At this point the average radiative force on the gas elements be- 
comes super-Eddington again, and so the gas at this radius expe- 
riences a net force outwards. A wind will therefore be accelerated 
from this radius. 

The forces acting on the material in this region are the pres- 
sure gradient outwards, and gravity inwards. Thus, the equation 
of motion of the gas (assuming a steady state wind, namely, that 
d/dt — > 0) is given by: 



dv 



GM 1 dP 



dr r'^ p dr ^ 

where v is the wind velocity. The continuity equation implies 

m = 47rr pv = const., 



(25) 



(26) 



where m is the mass loss rate. The equation of energy conservation 
is 



r, . ^ ■ I V^ / X GM\ 

L{r) = E -m[ h w{r) 



(27) 



where w — 5kT/2pmH + ^Prad/ p is the specific enthalpy and E 
is the total energy loss rate from the sonic point, given by 



Hj — J^sonic ~v rn 



— V, 

2 " 



+ w{Rsonic) ■ 



GM 



J^son 



(28) 



Because the wind is optically thick, the transfer of energy can be 
described by radiative diffusion. Note that although conditions for 
convection are still satisfied, we can neglect this energy transport 
since it is necessarily unimportant in the supersonic flow. Thus, the 
temperature gradient in the wind is 



dT 
dr 



3kLp 



16nacr^T^ ' 



(29) 



3.1.4 Photon Tired Winds 

The mass loss is predicted according to conditions at the sonic 
point. If the mass loss is too high and the potential well is too deep, 
the luminosity may be insufficient to push the wind to infinity. This 
gives rise to so called photon tired winds (Owocki et al. 2004). The 
mass loss for which the wind becomes photon tired is 



mti, 



LsRs 
MG 



(30) 
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The behavoir of photon tired winds was studied by van Marie et al. 
(2009). It was found that shocks form between infalling material 
and outflowing wind. This forms a layer of shocks in which there 
is a large kinetic flux, but without the associated mass flux. When 
photon tired winds arise, the mass loss from the top of the layer of 
shocks is reduced to less than the photon tiring limit, and we take 
it to be 

rh = 0.9mtiring. (31) 



3.1.5 Boundary Conditions 

The structure of the star is determined by the specific entropy at 
the center (which almost uniquely determines the mass) and by the 
central density, (see §1.1.1). The only free variable in a given struc- 
ture is therefore the luminosity. The value chosen for the luminos- 
ity should agree with the surface conditions on the radiation field, 
which should satisfy the blackbody radiation law: 



L 



ph 



47rrphCrTo'ff, 



where the Toff is the temperature at optical depth r 
optical depth in the wind is given by 

r ~ Kpr. 



(32) 
2/3. The 



(33) 



However, if the wind is sufficiently thick, the gas temperature can 
decrease enough for recombination of the ionized hydrogen to oc- 
cur before the photosphere is reached (i.e., at r > 1). This implies 
that the opacity decreases from its Thomson value to effectively 
zero, beyond which the gas is transparent. Thus, the photosphere in 
this case is going to be located at the recombination front. Since 
it is relatively sharp, we can take the classical Eddington result 
that r = 2/3 at this radius, and require eq. 32 to be valid there. 
The actual temperature of recombination depends on the density, 
but since the structure is not sensitive to the actual recombina- 
tion temperature, we take a nominal Tree ~ 4500K. As we will 
show in the results, the states of very low effective temperature, 
i.e., TcB < 4500K are very common in SMSs. 



3.2 Evolution of Non-Rotating SMSs 

Assuming that the central temperature never becomes high enough 
for nuclear burning to be important, and because in the mass range 
of 4 X 10'* to 10* Mq the thermal evolution timescale is longer than 
the hydrostatic timescale, the evolution of the SMS is described by 
a sequence of hydrostatic equilibrium states. The star begins as a 
large spherical cloud with little binding energy. As the star radiates 
its binding energy, it becomes denser and smaller. In addition, the 
wind from the star changes the total mass as the star evolves and 
the central entropy is no longer constant. For the evolution of the 
SMS, we need to Imow its total energy. The internal component is 
given by 



Eint = / 47rr udr, 



(34) 



where u is the internal energy density. The gravitational Newtonian 
energy is given by 



3.2. 1 End Stage of Non-Rotating SMS 

The star ends its hydrostatic equilibrium life when it reaches 
the critical radius, once the general relativistic corrections ulti- 
mately destabilize the star. This critical radius is given by (see 
Shapiro & Teukolsky 1983 for details) 



GM 



= 0.6295 



M„ \ 1/2 



\ M ) 



(36) 



AnrGmpdr. 



(35) 



RcritC^ 

Specifically, the criterion for instability is given in the form of a 
critical density (eq. 8) above which the destabilizing GR correc- 
tions overcome the stabilizing effects of the gas pressure. When 
this density is reached, the evolution of the star ends and a black 
hole is formed. 



3.3 Evolution of Rotating SMSs 

While the initial stages of rotating SMSs are the same as the non- 
rotating SMSs, the following steps change drastically. As men- 
tioned in §1.1.2, SMSs having enough angular momentum do not 
collapse when they reach the critical density. Instead, they keep 
evolving through a sequence of hydrostatic states, while radiating 
their binding energy and evaporating some of their mass. While 
contracting, the central temperature increases, and nuclear burning 
becomes important. At some point, the energy released by nuclear 
burning becomes large enough to support the SMSs from further 
contraction. 

We divide the evolution of rotating SMSs into 3 stages: 

• Stage la, Super-Eddington contraction: The evolution is gov- 
erned by the SMS contraction, and the radiation of its binding en- 
ergy. For non-rotating SMSs, this stage ends with GR collapse, 
which sets the duration of this stage. For rotating SMSs, the stage 
ends when the central temperature is high enough for hydrogen to 
be ignited. 

• Stage lb, Super-Eddington contraction with C buildup: Al- 
though the pp-chain and triple-a burning cannot supply a sufficient 
amount of energy to halt the contraction, the temperature can be 
high enough for these reactions to build up ^^C. This build up sets 
the duration of this stage. 

• Stage II, Super-Eddington CNO burning: The accumulated C 
slowly increases the fraction of energy produced by the CNO- cycle. 
Once it becomes dominant, the contraction stops and energy pro- 
duction proceeds through CNO burning. Like the previous stage, 
this one is super-Eddington and therefore has a continuum driven 
wind. As a consequence, the evolutionary timescale is governed by 
mass loss from the atmosphere. As the mass of the SMS decreases, 
so does its luminosity. 

• Stage III, sub-Eddington, "normal" massive star: At this 
stage, the luminosity is sub-Eddington, and the SMS becomes a 
"normal" massive star. Here, the evolutionary timescale is governed 
by nuclear burning. This final stage is the longest. Since the SMSs 
at this stage are normal stars, they have already been extensively de- 
scribed in the literature (e.g.. Bond et al. 1984), and therefore will 
not be described here. 



4 NUMERICAL METHODS 

The problem we are required to solve can be divided into two parts. 
First, we need to solve the hydrostatic model of the SMS, with a 
given mass and energy. Then, we need to evolve these models to 
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describe the evolution of the SMSs from their formation to their 
collapse. Here we describe the main methods we used. 

4.1 Solving a Stellar Model 

The stellar structure is obtained by solving eqs. 14-17 for the con- 
vective region, eqs. 19 and 23-24 for the atmosphere and eqs. 25-27 
and 29 for the wind. These equations are solved by first guessing 
the total luminosity. The solution is then iterated using the shoot- 
ing method, until the outer boundary condition given in eq. 32 is 
fulfilled. 



4.2 Evolving the Stellar Models 

After we obtain a single snapshot of a SMS, we can proceed to find 
how it evolves with time. Each model solution predicts a luminosity 
L and a mass loss rate m. The stellar mass and total energy are then 
evolved through conservation of mass and energy: 

dM _ . 

~dr ^ ~^ 

—7— ^ Lnuc — -tv, (38) 

at 



(37) 



where Lnuc is the total nuclear energy generation rate, while L is 
the total luminosity. E is the sum of the internal and gravitational 
energies. 

Since pc and s are the real variables describing the models, 
and not M and E, the model evolution proceeds by choosing a 
small central density increment Ape = epc (specifically with e — 
0.001), and then solving for the new model. 



5 RESULTS 

We now proceed to describe the results of our simulations. We 
present the structure and evolution of the SMSs from their early 
phase to the time they reach the critical radius and collapse, if they 
are non-rotating, or until they become normal massive stars, if they 
do rotate. 

In particular, we perform a parameter study. Since the exact 
"porous" opacity law required is not adequately known, we depict 
solutions of the same mass, but different effective opacity laws. 
This allows us to study the influence of the atmospheric effective 
opacity on the stellar structure and evolution. 



5.1 Quasistatic Configuration 

The early phase (-Rsonic ~ 2000i?t:rit) of a SMS with Afo ~ 
5.5 X 10'' M0 is depicted in figs. 2-3. In fig. 2 we present the gas 
density and the temperature in both the convective zone and the 
atmosphere. In the latter, the density decreases by three orders of 
magnitude, enabling a much "lighter" wind. 

The wind velocity and temperature are described in fig. 3, 
from the sonic point to the photosphere. In the case depicted here, 
the wind is optically and geometrically thick (Rph ~ 1.57?sonic). 



5.2 Stage la: Temporal Evolution Before Collapse or Ignition 
of Nuclear Burning 

Figs. 4-5 plot the temporal evolution of the sonic and photospheric 
radii, the luminosity and the mass of a star with an initial mass of 



A'l = 5.5 X 10''A/q. This star loses about 13% of its initial mass 
during this stage (until ingnition of nuclear burning). As can be seen 
in fig. 4, the wind is optically thick, causing the effective tempera- 
ture to remain low. Thus, a negligible amount of ionizing radiation 
is emitted during this stage of evolution, unlike the standard lore 
for SMSs. 

The results for non-rotating SMSs (which only have this stage 
of evolution) are summarized in table I , for different initial masses . 



5.3 Ignition and Evolution of Rotating SMSs 

Figs. 6-7 show stages lb and II of a rotating SMS with an initial 
mass M = 5 X lO'''A/0. Specifically, fig. 6 depicts the mass, lumi- 
nosity (as Fs), central temperature (Tc) and effective temperature 
(Tph) of this SMS. Fig. 7 describes the sonic and photospheric radii 
as a function of time. The two stages lb and II discussed in §3.3 
are easily distinguished. Stages lb lasts ~ lO'^yr. Here the SMS 
radiates its binding energy, while contracting (as can be seen as a 
fast decrease of the sonic radius). The central temperature increases 
from ~ 10^ to > 10* and the luminosity is very high (> lOLEdd)- 
Stage II lasts for ~ lO^yr. Here the evolution is ruled by the wind 
mass loss from the atmosphere. The total mass and luminosity de- 
crease, but the central temperature and the sonic radius remain al- 
most constant. Stage III (not shown) is the last and longest. At this 
stage, the luminosity is sub-Eddington and the evolution is ruled 
by the nuclear reactions. If a wind is present, it will be a weak line 
driven wind due to the small amount of CNO that was produced. 



5.4 The effects of the atmospheric opacity law 

One of the uncertainties in the model, is the opacity law behav- 
ior of the porous atmosphere. The reasons it is not known well is 
because it depends on the nonlinear radiative hydrodynamic config- 
uration the atmosphere will reach, and unfortunately, there are still 
no numerical simulations or empirical data which can constrain the 
effective opacity law. It is for this reason that we parameterized 
the effective opacity using eq. 21, with one free parameter, B. We 
have examined the influence of different values of B on the stellar 
structure and evolution. The values for B taken were B = 0.75, 
B = 1.0 and B = 1.25. 

We find that there is no significant difference in the appearance 
of SMSs, under different opacity laws. However, there are some 
differences with respect to SMS evolution. 

Because the luminosity depends on the opacity, the lifetime 
somewhat depends on the opacity law. For example, the lifetime 
of stages la and lb of a 5.5 x 10*A/q SMS (i.e., the lifetime of 
the gravitational collapse until ignition of CNO nuclear burning), 
ranges between 2800 yr for B = 0.75 and 7000 yr for B = 1.25. 

However, there is one major characteristic which is sensitive 
to the opacity law, and it is the mass of the sub-Eddington star that 
will end the super-Eddington phase. The opacity law can be de- 
scribed by a Fcrit, which is the Eddington parameter below which 
the system becomes sub-Eddington. Using the Eddington quartic 
relation and the approximation that the SMS can be described by 
an n = 3 polytrope, it is possible to estimate the critical Mass for 
which the SMS will become sub-Eddington (see Owocki & Shaviv 
2012); 



M, 



A* ^ crit 



18.3M, 



01 



(1 — I'crit) 

where p, is the molecular weight. 



(39) 
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Figure 2. Gas density p(r) (solid) and temperature (daslied) vs. stel- 
lar radius at tlie convective zone and tlie atmospliere for a SMS with 
5.5 X 10* Mq, in its early pliase (i?sonic ~ 2000iJcrit ~ 6 X lO^^cm). 



Figure 3. The wind velocity (solid) and the wind temperature (dashed) vs. 
radius for the same object as fig. 2. 



Since there is little depletion of hydrogen during the SMS 
evolution, we can take the primordial value of X ~ 0.75. More- 
over, Shaviv (2001a) found two instabilities, one with Tcrit ~ 0.8, 
which is more general, and one with Fcrit ~ 0.5, which de- 
pends on the boundary conditions. Other instabilities, when mag- 
netic fields are present could operate at even lower values of Tcrit 
(see for example the numerical analysis of instabilities in luminous 
accretion disks by Turner et al. 2005). For the above value of X 
and range of Fcrit, we find Merit ranging between 150 and 1200 
Mq . Namely, the "normal" sub-Eddington star left once the super- 
Eddington phase is over is very sensitive to the critical Eddington 
parameter above which atmospheres become porous. 



6 SUMMARY & DISCUSSION 

We obtained numerical solutions for rotating and non-rotating 
SMSs. In the range of masses 5 x 10'* to 1 x IQ^ Mq. In this 
range, the hydrodynamical timescale is shorter than the thermal 
timescale, such that the star is always in hydrostatic equilibrium 
during its evolution. Above these masses, there is no equilibrium 
phase and the star collapses on a dynamical timescale. For SMSs 
with M > 1 x IQ'^Mq, the thermonuclear energy generation 
of the non-rotating SMSs is negligible, thus, in this range, the 
SMSs are hydrostatic objects which slowly release their binding 
energy, and shrink until they end their life through a collapse into 
a super massive black hole. The rotating SMSs will not collapse, 
but instead ignite nuclear burning and lose mass through a wind, 
until a sub-Eddington star is left. The mass of this star depends 
on the Eddington parameter Fcrit for which the system becomes 
sub-Eddington (see §5.4). 

Our analysis included two main parts: 

• We constructed hydrostatic models of SMSs. These included 
the full calculation of the SMS structure while considering that lu- 
minous atmospheres are unstable and become "porous" as the Ed- 



Log(t[yrsJ) 

Figure 4. The sonic radius Rsonic (solid), and the radius of the photosphere 
ijph (dashed), as function of time for a star with an initial mass of M = 
5.5 X 10'' Mq. 



dington luminosity is approached, and that continuum driven winds 
are accelerated by these objects. When relevant, it also included nu- 
clear reactions. 

• We then evolved the models for SMSs using the hydrostatic 
configurations as basic building blocks. The evolution allows us to 
calculate various general properties such as the total mass loss, the 
energy radiated in ionizing photons, or the total life time. 

The main results obtained are the following: 

• The polytropic solution is a very good approximation where 
convection is dominant, but a very poor approximation in the outer 
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Mnitial 

[Me] 


rcollapaeilifetime) 

[yr] 


M(T)/A/i„itai 


Tc{r) 

[K] 


Pc{r) 
gr cm~^ 


5.5 X 10"* 


— 


0.87 


2.1 X 10** 


6 


1.0 X lO'"^ 


— 


0.86 


2.2 X 10** 


4.9 


2.5 X 10^ 


8.6 X 10^ 


0.94 


1.0 X 10* 


0.3 


5x 10^ 


3.2 X 10^ 


0.97 


5.2 X 10'^ 


2.7 X 10-2 


1.5 X 10" 


60 


0.98 


1.7 X 10^ 


5 X 10-* 


5.0 X lO^ 


10 


> 0.99 


5.0 X 10** 


8 X 10-"^ 



1.0 X lO'^ 



< 10 



1.00 



2.5 X 10^ 7.5 X 10"^ 



Table 1. Summary of the pre-coUapse stage of non-rotating SMSs . Tabulated are the various SMS characteristics as a function of the initial mass. A few 
points to note: (A) The two lower masses represented (5.5 X 10*Mq and 1.0 X 10^ Mq) will not collapse even if there is no rotation, because at some 
point of their evolution the thennonuclear energy generation becomes important, supporting the star against additional contraction. Thus, these objects are not 
SMSs, but "Very Massive Stars" instead (e.g., see Bond et al. 1984) (B) The SMSs always have a large enough photosphere, hence, a low effective temperature 
(determined by recombination in the wind). As a consequence, SMSs will not contribute a significant amount of ionizing radiation. (C) The mass shed by the 
objects is small relative to the total mass of the star. Thus, as long as thermonuclear reactions are not ignited, the total wind mass loss does not significantly 
affect the evolution. (D) Masses larger than the above range may collapse to become a black hole, but they do not forni an intermediate hydrostatic object 
because of the short gravitational collapse time scale. 





LogfflyrsJ) 



Log{t[yrsJ) 



Figure 5. The evolution of the mass (solid), luminosity at the base of the 
wind Ls/iEdd (dashed), central temperature Tc (dotted) and effective tem- 
perature Tpii (long dashed) as a function of time, for a star with an initial 
mass of M = 5.5 X IO^Mq. The SMS loses about 13% of its mass un- 
til ingnition of hydrogen. Note that Tpjj stays low during this stage, so no 
ionizing radiation is emitted. 



Figure 6. The evolution (stages Ib-II) of a rotating star with an initial mass 
Af = 5 X 10''' A/0. Plotted are the mass (solid), luminosity at the base 
of the wind (Ls/i'Bdd = Ts, dashed), central temperature (dotted) and 
effective temperature (long dashed). The effective temperature remains low 
during all stages of evolution, this SMS emitts very small amount of ioniz- 
ing radiation. 



parts of the star. We also note that some formation scenarios can 
produce a significantly non-polytropic structure (Begelman 2010), 
but these scenarios were not discussed here. 

• The solution describes objects which radiate above their Ed- 
dington limit, in contrast to the standard solution where the lumi- 
nosity is just below the Eddington limit. The luminosities obtained 
vary between ~ 2jCEdd to ~ 50£Edd at the base of the wind (the 
larger the mass the higher the luminosity), but the actual luminosi- 
ties emitted from the photosphere vary between ~ 10^'^/lEdd to 



~ 2i3Edd (the larger the mass, the lower is the photospheric lumi- 
nosity). The rest of the energy is used to drive a wind. 

• Because the binding energy at the onset of the GR collapse 
into a black hole does not depend on the outer layers of the star, 
higher luminosities imply shorter lifetimes in non-rotating stars. 

• The mass loss driven by the wind in the pre-collapse phase, is 
typically between 5% to 20% of the total mass, where less massive 
SMSs have larger mass loss. This immediately implies mass loss 
will be crucial in stars less massive than the least massive SMSs. 
These stars, which have prolonged lifetimes (due to nuclear reac- 
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the gas and the hmited efficiency of convection at its upper bound- 
ary, both effects are not due to the super-Eddington nature of the 
solution. 

Another intrinsic characteristic of the present solutions is the 
high mass loss rate. This has several effects. First, the high mass 
loss rates imply that large amounts of energy are required to push 
the gas out of the gravitational potential well. As a result, the lumi- 
nosity left above the wind can be significantly reduced, and even 
make the intrinsically super-Eddington SMSs appear to be sub- 
Eddington. 

Second, it implies that the photosphere resides at large radii, 
inside the wind. As a consequence, the photospheric temperature is 
lower by an order of magnitude relative to the standard SMS solu- 
tion. As a consequence, only a negligible amount of ionizing radia- 
tion is emitted. This is in markable contrast with the standard SMS 
models, and it implies that it is unlikely that SMSs have played a 
major role in reionizating the early universe. 



Log{t[yrsJ) 

Figure 7. Sonic radius rsonic (solid), and radius of the photosphere 
(dashed) of the same SMS as fig. 6, during stages lb and II of the evolu- 
tion. In stage III (not shown), the luminosity is sub-Eddington, so there is 
no continuum driven wind and the photosphere is located at the top of the 
atmosphere. 



tions), will lose significant amounts of mass, making mass loss the 
dominant factor in their evolution. 

• The effective temperature of the SMSs remains low during 
their evolution (i.e., the photosphere is determined by recombina- 
tion at the wind, taking place at ~ 4500K) due to the large photo- 
spheric radius. Hence, the amount of ionizing radiation emitted by 
these objects is negligible. 

• Rotating SMSs will lose most of their mass in the later stage 
(II, see §3.3), through their winds. The mass loss continues until a 
sub-Eddington massive star is left. 



6.1 Differences from the Standard Solution 

As mentioned before, the solution for SMSs obtained here is no- 
tably different from the standard lore of SMSs. It is thus worth 
while to emphasize these differences and their causes. 

In SMSs in general, the radiation to gas pressure ratio is very 
large. If the opacity is not modified and instead remains at its micro- 
scopic value of Thomson scattering, it is straightforward to show 
that the luminosity should be that of Eddington minus a small rel- 
ative coiTection of order the gas to total pressure ratio. Hence, the 
luminosity in the classical SMS solutions is taken to be the Edding- 
ton luminosity. 

However, this does not consider the rise of instabilities when 
the radiation pressure is dominant. Once taken into consideration, 
the effective opacity is reduced and the effective Eddington lumi- 
nosity is increased. This gives rise to solutions with L > £Edd, in 
which the outer regions of the hydrostatic part of the star is porous. 

Since the structure of the convective zone does not depend on 
the luminosity, this inner region (which comprises most of the stel- 
lar mass), does not change in the super Eddington solution. Never- 
theless, this region does not fit exactly the polytropic solution with 
a polytropic index of n = 3. This deviation from the polytropic so- 
lution is caused by two reasons: the finite pressure contribution by 
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